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and modification and its allied resources have progressively become a severe problem 
J prsscoes pulling the worldwide attention and now it rests at the central point of the 

conservation of the environment and sustainability. The present research aimed to 
examine the land-use changes and their impact on agricultural production using remote 
sensing and GIS techniques over the study area that comprised of Tehsil Shorkot, District 
Jhang, Punjab, Pakistan. Images were pre-processed by using the Arc GIS and ERDAS 
Imagine 15 software for stacking of the layers, sub-setting, and mosaicking of the satellite 
bands. After the pre-processing of the images, supervised image classification scheme was 
applied by employing a maximum likelihood algorithm to recognize the land-use changes 
which have been observed in the area under study. The area under water was occupied 9.6 
km* in 2010 that increased to 21.04 km* in 2015 and decreased to 19.4 km/*in 2020. Built-up 
land was 16.6 km? in 2010 that increased to 19.4 km’ in 2015 and 26.8 km?’ in 2020. The total 
area under vegetation was computed as 513.2 km*in 2010 that increased to 601.6km* in 2015 
and further increased to 717.7 km*in 2020. Forest land use showed decreasing trend as the 
covered area in 2010 was occupied 90.8 km* that decreased to 86.7 km*in 2015 and further 
decreased to 61.84 km’ in 2020. In 2010, barren land use was occupied 528.54 km* that 
considerably decreased to 429.64 km*in 2015 further decreased to 333.1 km*in 2020. Barren 
land drastically decreased into watered, built-up, and vegetation land uses. The findings of 
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this study will be helpful for the future conservation of various land-use types, urban and 
regional planning, and an increase in agricultural production of various crops in the study 
area. 

Keywords: Landuse; GIS; RS; agricultural production; Shorkot. 


Introduction 

Land modification and its allied resources have progressively become a severe 
problem presently pulling worldwide attention and now it rests at the central point of the 
conservation of the environment and sustainability [1]. It is perhaps due to the growing 
landuse manifestation and change consequences [2, 3]. The resulting impacts have been 
more emphasized by the enhanced land-use changes as a result of increased growth in 
population and the associated development in the socio-economic aspects from the different 
perspectives [4, 5, 6]. It is a general perception that nations attain economic prosperity 
having an increased population growth and affect the patterns of the consumption of many 
people, thus the increased needs by the well-off urban residents for the meat, nutrients, 
food, and dairy products have left untold negative impacts on the natural resources and also 
on the land and agricultural production [7]. A large number of models are in use for 
predicting the land-use changes for example CLUE [8], cellular automata (CA) [9], and 
Markov chain [10]. Land-use change concerning agricultural land is not only a role of the 
conditions [11] related to biophysical that makes a parcel of the land viable for agricultural 
production of various types [12]. Planning for urban land-use has especially attained much 
consideration in many studies and less attention has been paid towards the planning of 
acricultural land use [13]. Agricultural land is highly at the risk of conversion to non- 
agricultural undesired uses (housing colonies, roads, etc.) in many cities of the developed 
world as it has been witnessed in Bahawalpur city of Pakistan [14,15,16,17,18]. So, it is highly 
desired for efficient and more effective planning for the use of the agricultural lands in 
future. The basic motivation behind conducting this research was to formulate a network 
that permits the exploration of the substitute ways for the development of agricultural land 
use and its production in the future. Land-use changes play an important role in maintaining 
a global ecosystem balance [19], the human variables for example government policies and 
the environment largely impacts the land-use changes [20, 21]. A better understanding of the 
land-use changes can be acquired from the past practices, present patterns of the land use, 
and the future land use trajectory [22, 23, 24]. 

Based on the pixel, the change detection technique has been categorized into three 
stages [25], its features and the levels of the objects for processing of images and offers a 
classification of land use by using various techniques to detect the change [26,27,28]. It has 
been revealed from the past studies that the land-use changes pose a prospective to possess 
an important natural impact that changes the agricultural landuse [29]. Various studies have 
also endorsed the effects of the environment on land-use changes [30, 31, 32, 33]. Remote 
sensing (RS) and Geographic information system (GIS) have been very beneficial and 
important for evaluating and examining land-use changes [34, 35]. Remote sensing based on 
the satellite has the capability to offer synoptic land-use changes for a particular location and 
time period [36, 37, 38, 39]. Since the past few years, the studies have been largely dependent 
upon spatial remote sensing images data for mapping the individual species of the vegetation 


[40, 41, 42, 43, 44, 45]. 
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The significant objective of the current research was to examine the agricultural 
production impacted by land-use changes for the years 2010, 2015 and 2020 in the study 
area, with major classifications were incorporated te. water, built-up, vegetation, forest and 
barren land. Besides, the rice, wheat and sugarcane crops were incorporated to know the 
impact of temporal changes from 2010 to 2020. It was revealed that the water, forest and 
barren landuse area was decreased; in contrast, the built-up area was increased. 

Material and Methods. 
Investigation site. 

The word of the Shorkot has been derived from the word Shor, in Urdu, the 
sionificance of the word lies with the water logging and saltiness. This is the capital city of 
the Tehsil Shorkot, District Jhang, Punjab, Pakistan. Shorkot is situated between the 30° 30! 
N to 72° 24' E. This city is located on the Jhang-Multan road at an approximate distance of 
56 km from the district Jhang. Toba Tek Singh is situated in the NE with an approximate 
distance of 35 km. The elevation of the Shorkot from the sea level is about 131 meters 
(433ft). Shorkot occupied the status of Town Council before the year 2001. But later on, the 
city enjoyed the status of the Tehsil Municipal Administration (IMA). The area of tehsil 
Shorkot is approximately 1,158 square kilometers with an overall population of 670,255 
according to the 2017 census. Wheat, sugar cane and rice are the major crops sown in this 
area. The map of the study area is shown in figure 1. 


PAKISTAN 
TEHSIL SHORKOT 


PUNJAB 


DISTRICT 
JHANG 


Figure 1. Study Site 
Material and method. 

The landsat satellite images were acquired for the years 2010, 2015 and 2020 as 
mentioned in Table 1 from the website of United States Geological Survey (www.usgs.gov) 
and www.glovis.usgs.gov) and classified it into major land use classes. The crops data was 
collected from Punjab Agriculture Department, TehsilShorkot, District Jhang Office for the 
years 2010 to 2020. 

Image processing and classification (Maximum Likelihood) 

Landsat/Sentinel-2 images comprising of different bands were passed through 
vatious processes e.g., layer stacking in software ERDAS imagine 15 (stacking is a technique 
to generate the multiband image from a distinct band) mosaicking (to combine the two 
discrete images) and sub setting (study area was extracted after stacking of the image) [46]. 
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These composites were utilized to run the maximum likelihood classification algorithm to 
classify the land use and to demarcate major changes in the study area. Training samples 
were chosen by delineating the polygons for each land-use type by using the Arc GIS 10.8 
[47, 48]. Spectral signatures along with ground validation through Differential GPS were 
performed to apply training samples in supervised classifications. A supervised classification 
method was used to classify the land use images, the characters of these remote sensing 
images have been provided in Table 1. Pairs of geostatistical modeling and multivariate 
techniques were used to predict and examine the Spatio-temporal tendency of land-use 
change impact on agricultural production. 


Table 1. Characteristics of acquired remote sensing images 


St No Satellite Pixel-Size Year 
1 Landsat 4-5 30m 2010 
2 Landsat 8 15m (Pan-sharpened) 2018 
3 Sentinel-2 10m 2020 


Classification Accuracy Assessment 

The overall accuracy for a particular image or a classified map can be obtained by 
calculating the classification agreement of the Cohen’s Kappa coefficient [49]. Kappa 
statistics was used to assess and evaluate the overall classification accuracy. Kappa index was 
prepared to evaluate the classified maps to find out the accuracy of the results for each 


image. It is used to tolerate the change effects in the accuracy, its value should not fall below 


0.85. 


os Number of Correctly Classified Pixels in a class 
AS See. 
y Total Number of Pixels ina Class 


Number of Correctly Classified Pixels in a class 
Producer Accura¢gy = —#W7——_————- 
Total Number of Pixels in all Classes 


Total Number of All Correctly Classified Pixels 
Overall Accuracy = “wold wher of Be dlicaaa 
Classification accuracy assessment 

For classification accuracy, 100 reference points were generated. Each point had a 
pixel value, which was automatically identified by the software. Randomly generated points 
were then identified, and the corresponding class was assigned by the user manually. The 
error matrix and kappa statistics for the classified images were generated. The error matrix 
indicates the accuracy of classification [50] where rows represent the classes resulting from 
the classified image, whereas the columns represent the classes identified by the user from 
the reference values. The diagonal cells of the error matrix indicate the total number of 
correctly identified pixels for each class of the reference and classified data. The off-diagonal 
cells represent incorrectly identified pixels, which indicate the error between reference data 
and classified data. There are two types of errors, namely omission and commission error, 
are occurred during the classification process. 

Errors of commission occurred when a classification process assigns pixels to a 
specific class that doesn't belong to it. The number of pixels that were mistakenly assigned to 
a class was found in column cells of the class above and below the main diagonal. The 
Producer’s accuracy also described the number of errors of commission. For every class, 
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errors of omission occurred when pixels that belong to one class, are included in other 
classes. In the confusion matrix, the number of omitted pixels was found in the row cells to 
the left and the right from the main diagonal. The user's accuracy is another indicator 
characterizing the errors of omission. 

Results and discussion 


The land-use changes are very critical in the study area due to unintentional and 
unplanned built-up areas. Such changes in land use resulted in decreased forest and 
vanishing barren land and an increase in the built-up area. The change in study sites from 


2010, 2015, and 2020 are mapped in Figure 2. 
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Figure 2. Land-use Change Analysis for the years 2010, 2015, 2020. 

Table 2 reveals that the area under water was occupied 9.6 km*in 2010 that increased 
to 21.04 km’in 2015and decreased to 19.4 km*in 2020. Built-up land was 16.6 km? in 2010 
that increased to 19.4 km*in 2015 and 26.8 km? in 2020. The total area under vegetation was 
computed as 513.2 km? in 2010 that increased to 601.6km* in 2015 and further increased to 
717.7 km*in 2020. Forest land use showed decreasing trend as the covered area in 2010 was 
occupied 90.8 km* that decreased to 86.7 km* in 2015 and further decreased to 61.84 km/ in 
2020. In 2010, barren land use was occupied 528.54 km* that considerably decreased to 
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429.64 km*in 2015 further decreased to 333.1 km*in 2020. Barren land drastically decreased 
into watered, built-up, and vegetation land uses. The findings of this study will be helpful for 
the future conservation of various land-use types, urban and regional planning, and an 
increase in agricultural production of various crops in the study area. 


Table 2. Landuse change analysis of the study area for the Landsat images for the 
year 2010 to 2020 


Landuse Class Area 2010 km” Area 2015 km” Area 2020 km?’ 
Water 9.6 21.4 19.4 
Built-up 16.6 19.4 26.8 
Vegetation 513.2 601.6 717.6 
Forest 90.8 86.7 61.84 
Fallow and Bare Land 528.54 429.64 5354) 
Total 1158.74 1158.74 1158.74 


Land Use Change Analysis from 2010 to 2020 
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Figure 3. Landuse change analysis of the study area for the year 2010 to 2020 

The analysis of the study area revealed that the land-use changes were the major 
driving factors behind influencing the agricultural production of the area. Figure 3 provides 
analysis of Landuse change in the study area for the year 2010 to 2020. The rapid reduction 
in the forest, vegetation, and barren land provides a clear signal for the increase in the built- 
up and vegetation area. Another important objective of the present was to assess the impact 
of the land-use changes on agricultural production. Results showed that agricultural 
production of wheat, sugarcane and rice in the year 2010 was 34, 575 and 33 mounds (1 
mound= 40 ke) per acre respectively. Agricultural production of wheat, sugarcane and rice in 
the year 2020 was declined to 27, 900 and 36 mounds per acre respectively. The production 
and cultivable area of wheat crop decreased while the production and cultivable area of rice 
and sugarcane increased from 2010 to 2020 along with its cultivable area (Table 3). 

Table 3. The production of wheat, sugarcane and rice crops in the years 2010 to 2020 


S. Year Wheat Cultiva Sugar Cultiv RicePer Cultiva 
No Per Acte ble cane Per able Acre ble 
Productio Land Acte Land Productio Land 
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n (actes) Production (actes n (actes) 

(mounds) (mounds) ) (mounds) 
1 2010 34 184,580 oils 49,270 33 54,940 
Z 2011 34 185,980 600 50,240 34 56,350 
9 2012 23 182,350 600 49,550 34 56,900 
4 2013 34 179,990 650 51,600 35 57,200 
> 2014 34 175,290 700 52,250 ie 57,730 
6 2015 DA 175,550 700 52,503 34 59,195 
7 2016 oe 167,340 650 52,350 ys 60,250 
8 2017 of 168,250 800 53,450 ie, 59,540 
9 2018 one) 166,500 900 61,350 35 60,450 
10 2019 26 163,460 900 53,500 35 68,500 
11 = 2020 Ze 160,215 900 52,400 36 68,800 


Source: Agriculture Department, Tehsil Shorkot (2020) 
The confusion matrix were drawn to compute accuracy of classification that was 81%, 80% 
and 83% for the years 2010, 2015 and 2020, 
Table 4. Confusion matrix indicating the producer's accuracy, user accuracy, and overall 
accuracy 2010 LULC map of the study area. 


Reference Data 


, Built- Row sue 

Classified Data Soil Crop Vegetation Water Accuracy 
up Total (9/) 

Built up 29 6 0 0 0 35 82.86 
Soil 1 12 2. 0 0 iS: 80.00 
Crop 0 1 19 3 0 23 82.61 
Vegetation 0 1 3 14 0 18 77.78 
Water 0 1 1 0 7 9 77.78 
Column Total 30 21 25 17 fi 
Producer's 
Rea) 96.67 57.14 76.00 82.35 100.00 
Overall 
Classification 81 


Accuracy (%) 
Table 5. Confusion matrix indicating the producer's accuracy, user accuracy, and overall 


accuracy 2015 LULC map of the study area. 


Reference Data 


User 

Classified Data Built-up Soil Crop Vegetation Water 7 me Accuracy 
otal (%/) 

Built up a 3 0 0 30 90.00 
Soil 2 t 1 0 0 10 70.00 
Crop 0 Z 18 . 0 25 1200 
Vegetation 0 Q 4 15 1 20 75.00 
Water 0 1 1 0 13 15 86.67 
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Column Total 29 13 24 20 14 
— s Accuracy 9,149 53.85 75.00 75,00 92.86 
Overall Classification 20 


Accuracy (%) 
Table 6. Confusion matrix indicating the producer's accuracy, user accuracy, and overall 


accuracy 2020 LULC map of the study area. 


Reference Data 


User 


Classified Data Built-up Soil Crop ‘ee Water 7 we Accuracy 

otal 

(0) 

Built up 28 2 0 0 0 30 93.33 
Soil 3 17 0 0 0 20 85.00 
Crop 0 2 19 4 0 25 76.00 
Vegetation 0 0 4 10 1 15 66.67 
Water 0 1 0 0 9 10 90.00 
Column Total 31 22 23 14 10 
Producer's 
Acca Oo 90.32 77.27 82.61 71.43 90.00 
Overall 
Classification 83 


Accuracy (%) 
The scatterplots of results of classification are mapped in Figure 4. Confusion matrix 
indicating the producer's accuracy, user accuracy, and overall accuracy for LULC map of the 
study atea is presented in Table 3, 4, 5 and 6. 


Figure 4. Scatter plots. 
Conclusion. 

The study showcased that there was a rapid and unprecedented change in the barren, 
forest and water land-use class to built-up area and vegetation land-use. The decrease in the 
forest and barren land use classes were one of the main reasons for the reduction in 
agricultural production and increase in built-up area. The findings showed that historic 
conversion took place in the study area due to a decrease in the forest and barren land due to 
an accelerated increase in the density of population, development in infrastructure and urban 
advancements. Moreover, there had been a continuous change in barren and forest land-use 
class to built-up and vegetation land-use class. In order to increase the agricultural 
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production and economy of the country, a balanced check on the land-use changes is highly 
crucial. Hence, the present research will help in improving the capability of the government 
to suggest the various remedial measures for decreasing the impact of the land-use change 
conversion from agricultural/ forest land use to built-up area for sustainable management of 
land use. It will be also helpful in land use planning and in estimation of agricultural 
production in the future. 
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